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ABSTRACT 

This paper describes how active galactic nuclei can heat galaxy-cluster plas- 
mas by driving convection in the intracluster medium. A model is proposed in 
which a central supermassive black hole accretes intracluster plasma at the Bondi 
rate Meondi and powers a radio source. The central radio source produces cos- 
mic rays, which mix into the thermal plasma. The cosmic-ray luminosity L CT is 
oc MBondiC 2 . The cosmic-ray pressure gradient drives convection, which causes 
plasma heating. It is assumed that plasma heating balances radiative cooling. 
The plasma heating rate is self-regulating because Meondi is a decreasing func- 
tion of the specific entropy near the cluster center sq; if heating exceeds cooling, 
then s increases, M Bon di and L cr decrease, and the convective heating rate is 
reduced. This paper focuses on the role of intracluster magnetic fields, which 
affect convection by causing heat and cosmic rays to diffuse primarily along 
magnetic field lines. A new stability criterion is derived for convection in a 
thermal-plasma/cosmic-ray fluid, and equations for the average fluid properties 
in a convective cluster are obtained with the use of a nonlocal two-fluid mixing 
length theory. Numerical solutions of the model equations compare reasonably 
well with observations without requiring fine tuning of the model parameters. 

Subject headings: cooling flows — galaxies:clusters:general — galaxies: active - 
convection — magnetic fields — turbulence 



1. Introduction 

In many clusters of galaxies the radiative cooling time of intracluster plasma within 
the central 50 kpc is shorter than a cluster's age. In the cooling-flow model plasma heat- 
ing is neglected and the rate at which plasma cools to low temperatures (Mqf) is large, 
exceeding 10 3 M Q yr _1 in some clusters (Fabian 1994). Recent high-spectral- resolution x-ray 
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observations, however, show that the actual rate at which plasma cools to low tempera- 
tures is < 0.1M CF (Peterson et al 2001, 2003; Bohringer et al 2001; Tamura et al 2001; 
Molendi & Pizzolato 2001). This discrepancy, some times referred to as the "cooling-flow 
problem," suggests that plasma heating approximately balances radiative cooling in cluster 
cores. Plasma heating is also relevant to other unsolved problems in the study of galaxy 
clusters. Numerical simulations of cluster formation show that some heat source in addition 
to supernovae is needed to explain observed temperature and density profiles, star-formation 
histories, and mass-temperature and luminosity-temperature scaling relations (Suginohara & 
Ostriker 1998; Lewis et al 2000; Voit et al 2002; Tornatore et al 2003; Nagai & Kravtsov 2004; 
Borgani et al 2004; Dolag et al 2004). 

This paper focuses on plasma heating by active galactic nuclei (AGN) at the centers of 
clusters. The importance of AGN heating is suggested by the observation that almost all 
clusters with strongly cooling cores possess active central radio sources (Eilek 2004). Several 
mechanisms have been considered for transferring AGN power to intracluster plasma, includ- 
ing magnetohydrodynamic (MHD) wave-mediated plasma heating by cosmic rays (Bohringer 
& Morfill 1988; Rosner & Tucker 1989; Loewenstein, Zweibel, & Begelman 1991), Compton 
heating (Binney & Tabor 1995; Ciotti & Ostriker 1997, 2001; Ciotti, Ostriker, & Pelle- 
grini 2004), shocks (Tabor & Binney 1993, Binney & Tabor 1995), and cosmic-ray bubbles 
produced by the central AGN (Churazov et al 2001, 2002; Reynolds 2002; Bruggen 2003; 
Reynolds et al 2005). Such bubbles can heat intracluster plasma by generating turbulence 
(Loewenstein & Fabian 1990, Churazov et al 2004) and sound waves (Fabian et al 2003; 
Ruszkowski, Bruggen, & Begelman 2004a,b) and by doing pdV work (Begelman 2001, 2002; 
Ruszkowski & Begelman 2002; Hoeft & Bruggen 2004). Although considerable progress has 
been made, the coupling between AGN and cluster plasmas is still not well understood. 

The present paper describes how AGN can heat intracluster plasmas by driving convec- 
tion. In section 2 physical arguments are used to obtain the convective stability criterion in 
magnetized low-density plasmas containing cosmic rays. Section 3 develops a steady-state 
model of convective intracluster plasmas. In this model, a central supermassive black hole 
accretes intracluster plasma at the Bondi rate Meondi, leading to cosmic-ray production by a 
central radio source. The cosmic-ray luminosity L CT is taken to be proportional to MeondiC 2 ) 
and it is assumed that the cosmic-rays mix into the thermal intracluster plasma. The cosmic- 
ray pressure gradient drives convection, which is treated with a two-fluid (thermal-plasma 
and cosmic-ray) mixing length theory. Convection heats the plasma in three ways, by mix- 
ing hot plasma from the outer regions of a cluster in towards the center of the cluster, by 
providing a vehicle for cosmic-ray pressure to do work on the thermal plasma, and through 
the viscous dissipation of turbulent motions. It is assumed that plasma heating balances 
radiative cooling. The heating rate is self-regulating because Meondi is a decreasing func- 
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tion of the specific entropy near the center of the cluster, s ; if heating exceeds cooling, 
then s rises, Meondi falls, and the convective heating rate drops (Nulsen 2004; Bohringer 
2004a). The model leads to a set of coupled ordinary differential equations that are solved 
numerically using a shooting method. 

In section 4 the model is compared to observations. For a single choice of the model's free 
parameters, the solutions to the model equations compare reasonably well with observations 
of the Virgo cluster, the Perseus cluster, and Abell 478. However, the model underestimates 
the plasma density in the central ~ 30 kpc. This discrepancy may arise because the model 
does not properly account for the observed incomplete mixing of cosmic rays into the thermal 
plasma in the central region, a point discussed further in section 4. 

The model differs from an earlier two-fluid convection model (Chandran 2004 — hereafter 
paper I) primarily by including the effects of anisotropic conduction and cosmic-ray diffusion 
on convection. Anisotropic transport allows smaller cosmic-ray pressures to drive convection, 
which enables the model of the present paper to provide a better fit to the data. (In paper I, 
a larger central cosmic-ray pressure leads to an even smaller thermal pressure and density 
in the central ~ 30 kpc.) 

In addition to the topics mentioned above, section 4 also describes the relation of the 
model to other turbulent heating models and addresses the consequences of the model for 
diffuse gamma-ray emission from clusters and the mixing of heavy elements in the intracluster 
medium. 



2. Convective stability in a low-density plasma containing cosmic rays and a 

weak magnetic field 

In this section, physical arguments are used to obtain the stability criterion for convec- 
tion in a magnetized plasma/cosmic-ray fluid. Figure 1 depicts a fluid parcel rising from 
an initial radial coordinate r to a final coordinate r\ = r + Ar. Only instabilities on 
length scales <C r are considered. The gravitational acceleration g is in the — r direction. 
The horizontal solid line is a field line in the equilibrium magnetic field, which is taken to 
be perpendicular to g. The dashed line depicts how the field line is displaced as the fluid 
parcel rises. It is assumed that the magnetic field is strong enough so that heat and cos- 
mic rays diffuse primarily along magnetic field lines (a requirement easily met in clusters) 
but sufficiently weak that the Lorentz force can be ignored. The latter assumption may 
be too restrictive for clusters, but is maintained throughout this paper for simplicity. It 
is assumed that equilibrium quantities depend only on r. The equilibrium temperature, 
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thermal-particle density, and cosmic ray pressure at r and r\ are denoted (T , no,p cr ,o) and 
(Ti, ni,p crj i), respectively. Initially, the fluid properties in the parcel equal the equilibrium 
values at r . After the parcel is displaced to r 1; the fluid properties within the parcel are 
denoted (T f , n f ,p cr f). Requiring that the total pressure within the displaced parcel equal the 
total pressure in the surrounding medium gives 

n f k B T t + p CTif = n x k B T x + p CIjl . (1) 

If the frequency u of the fluid's linear modes were constrained to be either real or imaginary, 
then near marginal stability uj would approach and the time scale for the parcel's displace- 
ment to grow or oscillate would be very long. In this case there would be plenty of time for 
thermal conduction and cosmic-ray diffusion to smooth out T and p cr along the displaced 
magnetic field lines, giving Tf = T and p cr j = p cr ,o- Thus, if uj 2 were constrained to be real, 
then near marginal stability equation (1) could be rewritten 

n f k B T + p crfi = n x k B (t + Ar^j + Pcr,o + Ar ~^T' ( 2 ) 

where 7\ and p CXt \ have been expanded in a Taylor series about r = ro, and terms that are 
quadratic or higher order in Ar have been discarded. Equivalently, 

(n f - ni )k B T = Ar {n.k^ + ^ . (3) 

If ri{ — rii > 0, the displaced parcel is heavier than its surroundings and sinks back down, 
and the fluid is stable to convection. Conversely, if rif — n\ < 0, the parcel is buoyant and 
the fluid is unstable. (It is assumed that the mean molecular weight is a constant.) Thus, if 
uj 2 were real, the convective stability criterion would be 

dT dp cv 

nk B — + — — > 0. (4) 
ar ar 

[The subscripts can be dropped in equation (4) since Ar is taken to be small.] In a future 
publication a linear stability analysis will be used to show that equation (4) is in fact the 
necessary and sufficient condition for convective stability in the plasma/cosmic-ray fluid, 
even though uj 2 is not in general real. 1 In the absence of cosmic rays, equation (4) reduces 
to the stability criterion dT/dr > derived by Balbus (2001). [In stellar interiors, heat is 
transported primarily by photons, the conductivity is isotropic, and Balbus's criterion does 
not apply (Balbus 2001).] 



1 Equation (4) is the 7 — > 1 and 7 cr — ► limit of equation (28) of paper I, essentially because T and p c 
can be treated as constant in the displaced parcel near marginal stability. 
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In section 3, the above physical arguments are extended in an approximate way to 
galaxy clusters, in which the magnetic field lines are tangled. The same stability criterion 
is recovered, essentially because conduction and diffusion still act to resist changes in the 
temperature and cosmic-ray pressure of displaced fluid elements. Throughout most of a 
cluster, dT/dr is much less than T jr. On the other hand, for cosmic rays diffusing outwards 
from a central point source, dp cr /dr ~ —p C r/r. When dT/dr <C T/r and dp cr /dr ~ —p CT /r, 
equation (4) implies that cosmic rays can drive convection even when p cr <C p. Centrally 
produced cosmic rays in clusters are thus highly destabilizing to convection. 




Fig. 1. — A rising fluid parcel and one of the magnetic field lines threading the parcel. 



3. The AGN-driven convection model 

The intracluster medium is described using a standard set of two-fluid equations for 
cosmic rays and thermal plasma (Drury & Volk 1981, Jones & Kang 1990), modified to 
include thermal conduction, viscous dissipation, and radiative cooling: 

| = -„V.v, ,5) 

^ = _ 7pV • v + ( 7 - l)[# diss + V • (« • VT) - R], (7) 

and 

^ = -7crPcrV • V + V • (D • V Pcr ) + ( 7 cr - l)^inj, (8) 
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where 

d 9 „ 

p is the plasma density, v is the bulk velocity of the two-fluid mixture, p and p cr are the plasma 
and cosmic-ray pressures, T is the plasma temperature, $ is the gravitational potential, II 
is the viscous stress tensor, 7 and 7 cr are the plasma and cosmic-ray adiabatic indices (which 
are treated as constants), H diss is the rate of viscous heating, k is the thermal conductivity 
tensor, R is the radiative cooling rate, E in j is the rate of injection of cosmic-ray energy per 
unit volume by the central radio source, and D is an effective momentum-averaged cosmic-ray 
diffusion tensor. 2 

Cluster magnetic fields are easily strong enough to cause cosmic rays and heat to diffuse 
primarily along field lines, so that 

k ~ K\\bb, (10) 

and 

D~D{bb, (11) 

where b is the magnetic field unit vector, and k\\ and D\\ are the parallel conductivity and 
diffusivity. Cluster fields may also be strong enough that the Lorentz force and resistive 
dissipation are important [see, e.g., Pen, Matzner, & Wong (2003)], but for simplicity these 
terms are neglected in equations (6) and (7). The parallel conductivity is set equal to the 
Spitzer conductivity for a non-magnetized plasma, 



, , k-^T \ / 10~ 2 cm _3 \ / 37 \ cm 2 
K „ = k s = 9.2 x itfVfce f^tt — ~ hr*" — > (12) 



5/2 / IP' 2 cm' 3 \ / 37 \ — 2 
5 keVy V n o ) V lnA c 

where n e is the electron density, k-Q is the Boltzmann constant, and lnA c is the Coulomb 
logarithm. The value of D\\ is discussed below. 

The analytic fit of Tozzi & Norman (2001) is used to approximate the full cooling 
function for free-free and line emission: 



R = n{n c 



kuT\- 1 - 7 ( k-pT N 1 ' 



10" 22 ergs cm 3 s _i , (13) 



where n- x is the ion density and the numerical constants correspond to 30% solar metallicity. 
The mean molecular weight per electron \x e = p/n c mH is set equal to 1.18, and the ratio n;/n c 



2 Equations (7) and (8) are the same as the equations used in paper I, except that D and k are taken to 
be isotropic in paper I, and the definition of D in paper I corresponds to (j cr — times the value of D 
appearing in equation (8). 
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is set equal to 0.9. Radiative cooling of cosmic rays is ignored, which is reasonable if protons 
make the dominant contribution to the cosmic-ray pressure. Coulomb interactions between 
cosmic rays and thermal plasma are also neglected. 

The gravitational potential is taken to be the sum of a cluster potential and a central 
galaxy potential, 

$ = $ c + $ g , (14) 

with $ c corresponding to a fixed NFW dark matter density profile (Navarro, Frenk, & White 
1997), 

Pdm — — ( — : rz, {i-o) 

r[r + r s ) z 

where 

, 200 l_ 

° c 3 [ln(l + c) - c/(l + c)]' 1 J 

r s is the scale radius, c is the concentration parameter, and p crit = 3H 2 /8nG is the critical 
density at the redshift of the cluster (calculated assuming Q = 0.3, fl\,o = 0.7, and H = 
70 km s _1 Mpc _1 ). The central galaxy potential is set equal to 



1+ - 



. 2 

r 



(17) 



where v g and r g are constants, and t> g is the circular velocity for r ^> r g if $ c is ignored. 

Since the intracluster medium is turbulent, each fluid quantity is written as an average 
value plus a turbulent fluctuation: 

v = (v) + 6v, (18) 
P=(p} + Sp, (19) 

etc. The averaged quantities are taken to be spherically symmetric and independent of time. 
Equations (6) through (8) are then averaged. The viscous and inertial terms are neglected 
in the average of equation (6), so that 

S W> = -W^, (20) 

where 

Ptot=P + Pcv (21) 

The averages of equations (7) and (8) can be written 

y.F = -W+(H diss + V-(K-VT)-R), (22) 
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and 

V-F. = -Vy. + g^^ + <^), (23) 

where 

F = Fr = (24) 
7-1 1 J 

is the average internal energy flux, 

F cr = F cr r = *™zL, (25) 

7cr - 1 

W={pV-v), (26) 

and 

^cr = (PcrV-v). (27) 

The fluctuating quantities are treated as small, so (R) can be approximately obtained by 
replacing n e , n- u and T in equation (13) by their average values, and to a good approximation 

(p) = kB{p){T) . (28) 
jim H 

The mean molecular weight /i is set equal to 0.62. 

To obtain an average of the conduction and cosmic-ray diffusion terms, the disordered 
geometry of the magnetic field must be taken into account. Within the central 100 kpc of a 
cluster, the characteristic magnetic-field length scale Ib is probably in the range 1 — 10 kpc 
(Kronberg 1994, Taylor et al 2001, 2002, Vogt & Ensslin 2003, Ensslin & Vogt 2005), which 
is much smaller than the temperature gradient length scale Lt = (T) / (d{T)/dr). The field is 
also probably turbulent, with an inertial range of (shear- Alfven-like) magnetic fluctuations 
on scales ranging from ~ l B down to the proton gyroradius (Goldreich & Sridhar 1995, 
Quataert 1998, Narayan & Medvedev 2001). The chaotic trajectories of field lines lead to an 
isotropic effective conductivity kt for heat transport over distances ^> Ib, with kt reduced 
relative to ku (Rechester & Rosenbluth 1978, Chandran & Cowley 1998). The degree of 
reduction is not known with great accuracy, but approximate theoretical studies find that 
the reduction is by a factor of ~ 5 — 10 (Narayan & Medvedev 2001, Chandran & Maron 
2004, Maron, Chandran, & Blackman 2004). In this paper, it is assumed that 

kt = 9k\\, (29) 

where the constant 6 is a free parameter. The average of the conductive heating term is 
taken to satisfy the equation 



<v-(-vt)h-LA 



(30) 
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with T set equal to (T) in equation (12). Similarly, it is assumed that 

(V • (D • Vp cr )} = i A \r 2 D CI Up cr )] . (31) 
The value of D cr is taken to be 



D CT = y/Dl + vy, (32) 

where the constants Dq and are free parameters. The Vd term is loosely motivated by 
a simplified picture of cosmic-ray "self-confinement," in which cosmic rays are scattered by 
waves generated by the streaming of cosmic rays along field lines. If, contrary to fact, the 
field lines were purely radial, efficient self-confinement would limit the average radial velocity 
of the cosmic rays to the Alfven speed va, allowing the cosmic rays to travel a distance r in 
a time ~ r/v^- For constant v\, this scaling can be approximately recovered by taking the 
cosmic rays to diffuse isotropically with D cr oc r, the scaling that arises from equation (32) 
when Vdr D . This self-confinement scenario is too simplistic, since in clusters field 
lines are tangled, v a varies in space, and it is not known whether cosmic rays are primarily 
scattered by cosmic-ray-generated waves or by magnetohydrodynamic (MHD) turbulence 
excited by large-scale stirring of the intracluster plasma. It is not clear, however, how to 
improve upon equation (32). Self-confinement in the presence of tangled field lines is not 
well understood, and the standard theoretical treatment of scattering by MHD turbulence, 
which takes the fluctuations to have wave vectors directed along the background magnetic 
field, is known to be inaccurate (Bieber et al 1994, Chandran 2000, Yan & Lazarian 2004). 3 
A more definitive treatment must thus await further progress in our understanding of MHD 
turbulence and cosmic-ray transport. The value of D\\ is needed in the mixing length theory 
developed below. In analogy to equation (29), it is assumed that 

D\\ = (33) 



It is assumed that a black hole of mass Mbh at r = accretes intracluster plasma at 
the Bondi (1952) rate 

M = (34) 

corresponding to the average plasma parameters at r = 0, where c s is the adiabatic sound 
speed. The resulting values of M turn out to be a few M yr -1 , which is much smaller 



3 See Farmer & Goldrcich (2004) and Yan & Lazarian for a discussion of how MHD turbulence affects 
self- confinement, and Chandran ct al (1999), Chandran (2000b), and Malyshkin & Kulsrud (2001) for a 
discussion of the effects of magnetic mirrors. 
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than the mass inflow rate in the cooling-flow model (Fabian 1994). If one were to set 
(v r ) = —M/47rr 2 p with M of this order, the (v r ) terms in the average of equation (7) would 
be small compared to R for r > 1 kpc. Because of this, (v) is set to zero to simplify the 
model. 

The mass accretion rate inferred from equation (34) is used to determine the cosmic-ray 
luminosity through the equation 

L cr = r]Mc 2 , (35) 

where rj is a dimensionless efficiency factor. The spatial distribution of cosmic-ray injection 
into the ICM is not precisely known. Some clues are provided by radio observations, which 
show that cluster-center radio sources (CCRS) differ morphologically from radio sources in 
other environments. As discussed by Eilek (2004), roughly half of the CCRS in a sam- 
ple of 250 sources studied by Owen & Ledlow (1997) are "amorphous," or quasi-isotropic, 
presumably due to jet disruption by the comparatively high-pressure cluster-core plasma. 
With the exception of Hydra A, the CCRS in the Owen-Ledlow (1997) study are smaller 
than non-cluster-center sources, with most extending less than 50 kpc from the center of the 
host cluster (Eilek 2004). In this paper, the average cosmic-ray energy introduced into the 
intracluster medium per unit volume per unit time is taken to be 

(£inj> = S e~ r2 / r2 % (36) 

where the constant r s is a free parameter and the constant So is determined from the equation 
L cr = 47r f^drr 2 (E- m j) and equation (35). 

The F cr and W CI terms are related to F and W by the assumption that convection is 
subsonic. The perturbation to the total pressure is then neglected, so that 

Spcr = Sp. (37) 

Since (v) is also dropped, one has 
and 

W cr = -W. (39) 

The values of F and W are first obtained with the use of a local two-fluid mixing length 
theory that differs from the theory of paper I in that it accounts for heat and cosmic-ray 
transport along the magnetic field. Averages obtained in the local theory are denoted with 
the subscript L. These local quantities are then used as the basis for a nonlocal mixing 
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length theory. Average quantities in the nonlocal theory are weighted spatial averages of 
the corresponding quantities in the local theory, and are denoted with the subscript NL. 
The nonlocal quantities are used in the final model equations. The nonlocal theory is used 
because it is somewhat more realistic than the local theory and because it avoids the critical 
point discussed in paper I, which is related to the sensitivity of dF^/dr to the fluid profiles 
when the intracluster medium is near marginal stability. 

The local mixing length theory is derived as follows. Consider a fluid element in a 
convective region that starts at radial coordinate r ; at time t = with initial fluid properties 
Pi, Ti, p CTi i, and p to t,i equal to the average values at r = r x . Suppose that the fluid 
element undergoes a displacement d and is then mixed into the surrounding fluid, where 
|d-r| = I > 0, 

/ = ar (40) 

is the mixing length, and the positive constant a is a free parameter. For the moment, it 
is assumed that a < 1. Later, this assumption is relaxed — a well known inconsistency of 
local mixing length theory. The values of the fluid properties in the displaced fluid element 
just before it mixes into the surrounding medium are denoted pf, pf, T f , p cr f , and p to t,f- The 
rms value of 5v r is denoted u L . The approximate time for the fluid element to undergo its 
displacement d is 

At = — . (41) 
Integrating equation (5), one obtains 



Pi 

where 



* = e" r , (42) 



At 

dt V- v 

o 



(43) 

Lagr 



and 

i. At 



/•At 




/ dtf 


-/ 


/o 


Lagr J 



dtf(v(t),t) (44) 

is the (Lagrangian) time integral of the (arbitrary) function / evaluated at the location of 
the moving fluid element r(t). Since the fluctuating quantities and a are treated as small, 
|r| <C 1 and equation (42) can be written 

= -r. (45) 

Pi 
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It is assumed that the changes in p and p CT in the moving fluid element during the 
time At due to H diss , R, and E in j are small compared to the changes induced by the other 
terms in equations (7) and (8). Integrating equation (7) then yields 



Pf ~ Pi = - 7 



i-At 



dt pV • v 



Lagr 



At 



+ ( 7 -l) / dtV-(K-VT) 



(46) 



Lagr 



To lowest order in 5p/ (p) and a, 



At 



dtpV 



PiT. 



(47) 



Lagr 



When a fluid parcel of size / is displaced radially outwards, it remains magnetically connected 
to the same set of fluid elements to which it was initially connected (at least until it is 
mixed into the surrounding fluid, at which point it is assumed that the magnetic field in 
the fluid parcel is randomized). This is depicted schematically in figure 2. Because parallel 
conductivity acts to smooth out temperature variations along the magnetic field, the initial 
temperature of any fluid element in the parcel can be thought of as a weighted average 
of the temperature along the magnetic field out to some distance L». Since the parallel 
conductivity in clusters is considerable and a is treated as small, L\\ can be taken to be 
larger than /. If one ignores temperature changes occurring outside the fluid parcel, then 
parallel conductivity acts to restore the fluid element's temperature towards its initial value, 
and the last term in equation (46) can be approximated as 



( 7 -l 



pAt 

>/ ' 

JO 



dtV ■ (k • VT) 



Lagr 



£At( 7 - 1)^1 (Tf-Tj) 
I 2 



(48) 



where £ is a positive dimensionless constant, which is set equal to 1 in all of the numerical 
examples presented below. Equation (46) can thus be rewritten as 



Pi ~ Pi = ~ 7l>i - 



£At( 7 -l)K||(T f -Ti) 
P 



(49) 



Parallel cosmic-ray diffusion is treated in the same way as parallel thermal conduction. 
Integrating equation (8) thus yields 



Pcr,f Pcr,i 7crrPcr,i 



£At-D||(p CT[ f -p CT ,i) 

I 2 



(50) 



Since fluctuating quantities and a are treated as small, 

T f -7] Pi - Pi pi - pi 



Ti 



Pi 



Pi 



(51) 



Fig. 2. — Schematic diagram of a rising fluid parcel. The solid line is a magnetic field line passing through 
the parcel's initial location. The dashed line is an idealization of how the field line changes as a result 
of the fluid parcel's displacement. In the diagram, velocities far from the fluid parcel are ignored, just as 
temperature variations far from the fluid parcel are neglected to derive the estimate in equation (48). 



Equations (45), (49) and (51) combine to give 

*-»=(*- A)* fe^), (52) 

where 

£(7-lW 

di = : (53) 

lu L pi 

is roughly the ratio of At to the time for heat to diffuse a distance / along the magnetic 
field. When a\ <C 1 the thermal plasma expands adiabatically, and when a± >> 1 the thermal 
plasma expands isothermally. With the aid of equation (45), equation (50) can be rewritten 

y ' ' Pi(l+a 2 ) K J 

where 

a 2 = (55) 

is approximately the ratio of At to the time for cosmic rays to diffuse a distance / along 
the magnetic field. When a2 < 1, the cosmic rays expand adiabatically. When a 2 ^> 1 the 
cosmic ray pressure in the fluid element remains constant as the element is displaced, as in 
the linear Parker instability in the large- .Dm limit (Parker 1966, Shu 1974, Ryu et al 2003). 



Since it is assumed that the convection is subsonic, the total pressure in the fluid el- 
ement remains approximately the same as the average total pressure in the fluid element's 
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surroundings. Thus, to lowest order in a 

Ptot,f - Ptot,i = ^(Ptot)- (56) 

Since only the lowest-order terms have been kept, it is not necessary to specify in equa- 
tion (56) whether d(p to t)/dr is evaluated at n or n + /, so the argument of d(p tot )/dr is not 
explicitly written. Arguments of average quantities are omitted in similar instances below. 
Adding equations (52) and (54) and using equation (56), one obtains 

l-riPtot) = , (57) 

dr pi 

where 

A= (7 + ai)ft Tgg^ (58) 

1 + ai l + a 2 17 

The difference between the density inside the displaced fluid element and the average 
density of the surrounding medium is given by 

Ap = p f - (a + ^(p>) • (59) 
Upon using equation (57) for p f — p h one can rewrite equation (59) as 

A ' = '(^'>-»- < 60 > 

The fluid is stable if an outwardly displaced parcel is heavier than its surroundings. Since 
dA/du^ > and d{p tot )/dr < 0, the quantity in brackets in equation (60) is positive for any 
Ml > if it is positive as «l — * 0. The stability criterion is thus (where from here on angle 
brackets around average quantities are dropped) 

ldp^_l_dp >Q 
p dr p dr 

From a physical standpoint, the — > limit is appropriate for the stability boundary 
because the turbulent velocity is small near marginal stability. Since it is assumed that the 
mean molecular weight p is a constant, equation (61) can be rewritten 

dp cr , dT , . 

J T L +nk B —>0 62 
dr dr 



[as in equation (4)], where n = pj{pmu) is the number density of thermal particles. 
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If the stability criterion is satisfied, the local convective velocity ml is set equal to zero. 
Otherwise the fluid is convectively unstable and the value of u L is found by solving the 
fourth-order polynomial equation 



put 



lAp d$ 



16 dr 



(63) 



Equation (63) states that the mean radial kinetic energy of the fluid element is the mixing 
length times the buoyancy force on the fully displaced parcel times the numerical factor 
of 1/16 that is commonly used in one- fluid mixing length theory (Cox & Giuli 1968). Once u L 
is found, pf — pi and p { — p { are determined from equations (52) and (57). The difference 
between the pressure inside the displaced fluid element and the average pressure of the 
surrounding material is 

dp s 



Ap = p { - [pi + l 



dr 



(64) 



It is assumed that half of the fluid at any radius r consists of rising fluid elements and half 
consists of sinking fluid elements, each of which contributes comparably to the internal-energy 
flux. These fluid elements originated over a range of initial radii r\ in the interval (r — l,r+l), 
with pressure fluctuations whose magnitudes increase roughly linearly with \r— rj|. Averaging 
over rj, one obtains the estimate = c^gU-^Ap/ (7 — 1). The constant c avg is chosen to be 1/2 
to match standard treatments of one- fluid mixing length theory (Cox & Giuli 1968), 4 so that 
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F - Ml/ 


' p /7 + aA 


dptot 


dp 


L 2(7-1) 


A\l + aJ 


dr 


dr 


The typical value of V • v, denoted (div 


v), is given by 




(div v) = —— 
v 1 At 


/ d£V-v 
Jo 


Lagr 


u L T 
I ' 



(66) 



The value of Wl is estimated using the same factor of 1/2 employed in estimating F L , so 
that W L = Ap(div v)/2, or 



W, = — 



U L l 


P 




dptot 


dp 


dptot 


2A 


A 




dr 


dr 


dr 



(67) 



Note that in the limit (12 — > 00, A = p(j + ai)/(l + ai) and Ap = I dp CT /dr. This is because 
in this limit p CT) f = p CT ^. Thus, since ptot,f = Ptot(n + 1), Ap equals the change in the average 
cosmic-ray pressure of the surrounding fluid. 



4 In paper I, c avg was instead chosen to match values in the literature for the constant C\ defined in 
equation (85) below. 
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In the study of Ulrich (1976), the nonlocal heat flux is given by a weighted spatial 
average of the heat flux obtained from local mixing length theory. The same approach is 
adopted here, so that 

/oo 
d Zl F h ( Zl )ij F (z- zi), (68) 
-oo 



NL 

and 



oo 
oo 

Wi 



/oo 
dziW^z^wiz- Zl ), (69) 
-oo 



/oo 
dz l u h (z l )i) u (z - z^, (70) 
-oo 

where dz = dr/H p , and H p is the pressure scale height. In this paper, H p is replaced by r, 
so that 

^ ln fe)' (71) 

where r re f is an unimportant constant. Different forms for the kernel function ipp were 
considered by Ulrich (1976). Here, the following values are adopted: 

, t \ i ( \ f « _1 e^ /Q ifx>0 

M*) = M*) = { ifx < , (72) 

and 

a -} e -x/a w if x > 



Equations (68) through (73) are equivalent to the differential equations 

ar + F NL = F L , (74) 

ar — h«NL = «L, (75) 

dr 

and 

a w r + W NL = W L . (76) 
dr 

For r > r conv , where r conv is the largest radius at which the fluid is convective, F NL oc r~ x l a . 
When F and W are set equal to F NL and in equation (22), the term containing F NL 
is oc r -1-1 /" for r > r conv . To obtain the same scaling for the term containing Wnl in 
equation (22), the value of aw is determined from the equation 

a w l = a' 1 + 1. (77) 
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Equations (72) and (73) represent a one-sided average, in the sense that the nonlocal 
flux at radius r depends only on the fluid quantities at smaller radii. A more realistic two- 
sided average would correspond to second-order differential equations for F NL , w NL , and W^l 
(Travis & Matsushima 1973, Ulrich 1976). The one-sided average is adopted for convenience. 
[The reason is that the second-order differential equations have homogeneous solutions that 
grow with r, which make it more difficult to guess the boundary conditions at r = when 
the equations are solved numerically using the shooting method described below.] A more 
sophisticated nonlocal theory could be developed along different lines (see e.g. Ulrich 1976, 
Xiong 1991, Grossman, Narayan, & Arnett 1993), but is beyond the scope of this paper. 

The average of the viscous dissipation term is set equal to 

0.42p< ns 

-"diss — j , yio) 

where 

Urms = M NL (79) 

is the rms three-dimensional velocity corresponding to the rms radial velocity represented 
by «nl, and the constant 0.42 is taken from direct numerical simulations of compressible 
magnetohydrodynamic turbulence (Haugen, Brandenburg, & Dobler 2004). 5 [It turns out 
that in the solutions presented below, //diss is much smaller than the other turbulent heating 
terms (see footnote 4 of paper I for an explanation).] 

Setting F = F NL and W = Wnl, one can rewrite equations (20), (22), and (23) as 6 

dptot d<& 

^(^ = -^ L + ^ ss + ^(rVf)-^ (81) 



and 



Equations (74), (75), (76), (80), (81), and (82) form a system of four first-order equa- 
tions and two second-order equations for the six variables p, T, p CT , Fnl, «nl, and IVnl- 



5 The constant 0.42 is obtained by taking the mixing length I to correspond to ir/k p in the simulations of 
Haugen ct al (2004), where fc p is the wave number at which kE{k) peaks, and E(k) is the power spectrum 
of the turbulent velocity. 

6 See paper I for a discussion of the p cr — > limit and the need to include (v r ) to treat this limit properly. 



Eight boundary conditions are required to specify a solution. Two boundary conditions are 
obtained by imposing a density p ou tcr and temperature T outcr at radius r outcr . Two boundary 
conditions are obtained by taking dT/dr and dp cr /dr to vanish at the origin. Three addi- 
tional boundary conditions are obtained by requiring dF^/dr, dW^h/dr, and du^/dr to 
be finite at the origin. This leads to the condition at r = that Fnl = -Pl, W^nl = Wi,, 
and «nl = Uh- (Since dT/dr and dp„/dr vanish at r = 0, the cluster is marginally stable at 
r = 0, and F L , W L , and u L vanish at r = 0.) The eighth boundary condition is obtained by 
assuming that p CT — > as r — > oo. This condition is translated into a condition on p cr at r outcr 
as follows. The value of r oute r is chosen to be significantly greater than r s and much greater 
than D /vd, so that for r > r outer , E in j is negligible and D cr ~ t^r. 7 In addition, r outer is 
taken to lie outside r conv , the largest radius at which the intracluster medium is convectively 
unstable, so that F NL = F outcr (r /r^^)- 1 ^ and W NL = W^ out e r (r/r out er) -1-1/a for r > r outcr , 
where -F outcr and W outcr are the values of F NL and W^nl at r = r outcr . Solving equation (82) 
and requiring that p cr — > as r — > oo, one finds that for r > r outer 

dpcr _ X _ 2pcr 
dr Vdr 1+1 / a r 

where 

X = (2a - 1)( 7 - l)Foutc r rl Q cr + a( 7cr - 1)^^^. (84) 

Equation (83) applied at r = r outcr provides the eighth boundary condition. Numerical 
solutions are obtained using a shooting method. Values are guessed for p, T, and p CT at r = 
and the equations are integrated from to r out er- The guesses are then updated using 
Newton's method until the three boundary conditions at r outcr are met. 

Figure 3 compares three model solutions to x-ray observations of the Virgo cluster, 
Perseus cluster, and Abell 478. In all three solutions, the following parameters are used: 
M BH = 1O 9 M , r, = 0.01, 7 = 5/3, 7 cr = 4/3, a = 0.5, r s = 25 kpc, v d = 100 km/s, D = 
10 28 cm 2 /s, 9 = 0.2, v g = 400 km/s, r g — 1 kpc, and r oute r = 100 kpc. The only differences 
between the model solutions are the parameters describing the dark-matter density profiles, 
which are taken from the literature (see table 1), and the values of p outC r and T outcr , which 
are chosen to match the observations at r = r outcr in Perseus and A478 and at r = 61 kpc 
(the outermost data point) in the Virgo cluster. 

The deprojected temperature and density data for A 478 are provided by S. Allen 
(private communication; see Sun et al 2003) and assume Q = 0.3, Q^,o = 0.7, and H = 
70 km s -1 Mpc -1 . The deprojected temperature and density data for Perseus are provided 



7 The case Vd = requires a different approach and is not treated in this paper. 
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Fig. 3. — First, second, and third columns correspond to Virgo, Perseus, and Abell 478, respectively. 



by E. Churazov (private communication) and are based on the observations of Churazov 
et al (2003), but with H = 71 km s _1 Mpc _1 instead of the value H = 50 km s^Mpc -1 
used in the original article. The central peak in the Perseus temperature data is due to 
the hard power-law spectrum of the central active galaxy NGC 1275 (E. Churazov, private 
communication). The data for Virgo are from Matshushita et al (2002). 

Figure 4 shows the cosmic-ray pressure divided by the plasma pressure and the rms 
turbulent velocity for the three model solutions plotted in figure 3. For reference, in the 
model solution for the Virgo cluster, M = O.O81M yr _1 , L CT = 4.6 x 10 43 ergs/s, the 
radiative luminosity out to 100 kpc is 1.2 x 10 43 erg/s, and the radiative luminosity out to 
300 kpc is 2.5 x 10 43 erg/s. In the model solution for the Perseus cluster, M = 1.7M & yr _1 , 
L CT = 9.8 x 10 44 ergs/s, the radiative luminosity out to 100 kpc is 4.1 x 10 44 erg/s, and 
the radiative luminosity out to 300 kpc is 1.0 x 10 45 erg/s. In the model solution for the 
A 478, M = 3.8M yr^ 1 , L CI = 2.2 x 10 45 ergs/s, the radiative luminosity out to 100 kpc 
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Table 1: NFW parameters and redshift. 



Cluster r s (kpc) c z reference 



Virgo 560 2.77 (0.0040) 1 
Perseus 500 3.86 0.018 2 
A 478 610 3.88 0.088 3 



NOTE.— redshift (z) of Virgo is set equal to H d/c, with d = 17 Mpc and H = 70 km s _1 Mpc _1 . 
REFERENCES.— (1) McLaughlin (1999); (2) Ettori, Fabian, & White (1998) (3) Schmidt, Allen, & Fabian 
(2004). 

is 1.0 x 10 45 erg/s, and the radiative luminosity out to 300 kpc is 2.8 x 10 45 erg/s. 

4. Discussion 

Active galactic nuclei offer a promising solution to the cooling-flow problem on both 
observational and theoretical grounds: there is an active radio source at the center of almost 
every strongly cooling cluster (Eilek 2004), and AGN feedback in combination with ther- 
mal conduction can lead to globally stable equilibria (Rosner & Tucker 1989, Ruszkowski & 
Begelman 2002). The recent detection of moderately subsonic plasma motions in the Perseus 
cluster (Churazov et al 2004) suggests that turbulent heating is also important for the intr- 
acluster medium. The model presented in this paper connects AGN feedback to turbulent 
heating. It is assumed that radiative cooling is balanced by a combination of turbulent heat- 
ing and thermal conduction. The turbulence arises from convection, which is driven by the 
cosmic-ray pressure gradient. The cosmic rays are produced by a central radio source, and it 
is assumed that the cosmic rays mix into the thermal plasma. The model treats convection 
using a nonlocal two-fluid mixing length theory. In some ways the model is analogous to a 
convective star, with cosmic rays playing the role of photons. 

4.1. How successful is the model? 

The model is compared to observations in figure 3. For Perseus, the model temperature 
profile has a plateau-like feature at r ~ 20 — 30 kpc, which also appears in the observational 
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Fig. 4. — First, second, and third columns correspond to Virgo, Perseus, and Abell 478, respectively. 

data. This local flattening of the temperature profile makes the intracluster medium more 
convectively unstable, and is associated with large turbulent velocities. The model predicts 
maximum turbulent velocities of order 300 km/s in Perseus, similar in magnitude to the 
plasma velocities inferred by Churazov et al (2004). It should be emphasized that the three 
solutions shown in figure 3 are all based on a single choice of "round-number" values for the 
model's free parameters, as opposed to values that were fine-tuned to fit the data. Given the 
lack of fine tuning, the model is reasonably successful at matching the observations, although 
there are some important discrepancies. The model overestimates the temperature in A478 
at r ~ 200 — 300 kpc and underestimates the electron density at r ~ 10 — 30 kpc in each 
cluster. This density discrepancy persists for all values of the model's free parameters that 
have been investigated thus far and may be related to the assumption that the cosmic rays 
are well mixed into the thermal plasma, as discussed further in the next section. 
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4.2. Mixing of cosmic rays into the thermal plasma. 

The mixing- length theory developed in section 3 treats the turbulent fluctuations 5p, 
Sp, and 5p cr as small quantities. This entails the assumption that the cosmic rays are well 
mixed into the thermal plasma (since at fixed r there is little variation in p CT /p). However, 
roughly one-fourth of the clusters in the Chandra archive show clear evidence of x-ray cavities, 
which are typically ~ 10 kpc in radius at a projected distance of ~ 20 kpc from the cluster 
center (Birzan et al 2004). Many of these cavities are associated with enhanced synchrotron 
emission, suggesting that the cavities are cosmic-ray bubbles inflated by the central radio 
source. These cavities indicate that, at least within the central ~ 30 kpc of many clusters, 
the cosmic rays are not well mixed into the thermal plasma and 5p ~ (p), 5p CT ~ (p CT ), etc. 
It is likely that the cosmic rays are better mixed into the thermal plasma at larger r, owing 
to the combined action of cosmic-ray diffusion, instabilities of bubble/plasma interfaces, and 
turbulent mixing, but in this paper no attempt is made to quantify the degree of mixing as 
a function of radius. 

In low-p cr regions within the central 30 kpc, cosmic rays are less able to support the 
plasma against gravity than in the model of this paper, and \dp/dr\ has to be larger. Ac- 
counting for this properly would likely increase n e in parts of the central region, and possibly 
increase the plasma mass in the central region. Even if the plasma mass profile remained 
the same but the electrons were confined primarily to some fraction of the volume, the x- 
ray surface brightness would increase because R oc n 2 e . A larger surface brightness would 
in turn increase the observationally inferred value of the volume-averaged electron density. 
A proper treatment of incomplete mixing might thus resolve the discrepancy between the 
model density profile and the observations, but further work is needed to investigate this 
possibility. 

4.3. Relation to turbulent heating models 

Several studies have shown that for plausible values of the turbulent velocity and veloc- 
ity correlation length, heating from turbulent diffusion and viscous dissipation of turbulent 
motions can balance radiative cooling in clusters with a wide range of average tempera- 
tures (Kim & Narayan 2003, Voigt & Fabian 2004, Dennis & Chandran 2005). The main 
challenge for turbulent heating models is to explain how u Tms gets fine-tuned to the value re- 
quired for turbulent heating to balance radiative cooling. The AGN-driven convection model 
presented in this paper offers an explanation for this fine-tuning by linking the turbulence 
amplitude to the cosmic-ray luminosity L cr of the central radio source. The Bondi accretion 
rate Meondi is a decreasing function of the specific entropy of the central plasma, sq, and 
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L CT oc MeondiC 2 . Thus, if -u rms were so large that the total heating exceeded radiative cooling, 
then s would grow and M Bo ndi and L cr would decrease. This would in turn reduce p CT and 
thereby reduce M rms , since it is the cosmic-ray pressure gradient that drives the convection. 

4.4. Mixing of heavy elements 

Galaxy clusters with short central cooling times ( "cooling-core clusters" ) have centrally 
peaked heavy-element abundances (Fukazawa et al 1996, DeGrandi & Molendi 2001). While 
the iron abundance peaks sharply towards r = in such clusters, the abundance of oxygen is 
essentially flat in the central region. Since oxygen is effectively produced by type II SN but 
not by type I SN, the flat oxygen-abundance profile indicates that heavy element enrichment 
in cluster cores is dominated by type I SN. 8 For plausible models of the history of the SN la 
rate in the central CD galaxy, more than 5 Gyr is typically required for SN la to enrich the 
plasma in cooling-core clusters within the central ~ 50 kpc to the observed levels (Bohringer 
et al 2004b). This implies that turbulence can not transport metals out of the central 50 kpc 
on a time scale much shorter than 5 Gyr. At the same time, the abundance peak is spread 
out over a larger region than the stars of the cD galaxy. Thus, some mixing out of the central 
region is required (Bohringer et al 2004a). 

The eddy diffusivity in hydrodynamic turbulence is 

-Deddy = CiM rms /o, (85) 

where u rrns is the rms turbulent velocity and l is the dominant length scale (outer scale) 
of the turbulent velocity. Values in the literature for the dimensionless constant c\ range 
from 0.06 to 0.18 (see Dennis & Chandran 2005), although it should be noted that these 
values depend on the definition employed for l . [These values for c\ are <C 1 because in a 
3D random walk the diffusion coefficient is (Ax) 2 /6At instead of just (Ax) 2 / At, where Ax 
is the step size and At is the time between steps (Chandrasekhar 1943).] If Iq = ar, then the 
time scale for elements to be mixed out of the region interior to radius r, t mix = r 2 /Z) e dd y , is 

„ / r \ /200 km/s\ /0.1\ /0.5\ , . 

Some studies have argued that in MHD turbulence, magnetic tension reduces -D e dd y relative 



8 In the bulk of a cluster's volume, the iron-to-silicon ratio is typical of type II supernovae, and it is 
generally believed that type II SN dominate the heavy element enrichment (Finoguenov et al 2000, DeGrandi 
& Molendi 2001). 
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to the estimate in equation (85) (Vainshtein & Rosner 1991, Cattaneo 1994), but this claim 
has been recently challenged (Cho et al 2003). 

Equation (86) indicates that mixing times are comparable to enrichment times for plau- 
sible turbulence parameters, although a more careful analysis is needed to determine whether 
the AGN-driven convection model is consistent with enrichment-time constraints. Two fac- 
tors that would affect such an analysis are the following. First, the turbulent velocity peaks 
in the central region, dropping significantly outside of the central 50 kpc in the solutions 
plotted in figure 4. This allows for efficient mixing within the central 50-100 kpc, but reduces 
the rate at which metals are mixed out to larger radii. Second, the observed incomplete mix- 
ing of cosmic rays in the central ~ 30 kpc discussed above may limit the ability of rising 
cosmic-ray-dominated bubbles in the very central region to carry enriched plasma outwards 
(see, e.g., Briiggen et al 2002). 



4.5. Implications for diffuse 7-ray emission 

Cosmic-ray protons in the intracluster medium interact with thermal-plasma nucleons 
to produce pions. Neutral pions then decay into gamma rays. In this section, the results 
of Pfrommer & Ensslin (2004) are used to calculate the gamma-ray fluxes that would result 
from the model solutions for p cr (r) and n e (r) presented in section 3. The cosmic-ray-proton 
distribution function (number of cosmic-ray protons per unit volume per unit momentum), 
denoted /, is assumed to be 

/M-vM^P^. (87) 

where the constant a p is a free parameter satisfying 2 < a p < 3. At the risk of notational 
confusion, p in equation (87) is the cosmic-ray momentum, even though p was previously 
defined to be the plasma pressure. It is assumed that all of the cosmic-ray pressure arises from 
cosmic-ray protons. The value of n cr (r) corresponding to the model solutions of section 3 

POO 

is thus determined by equating the value of p CT {r) in the model solutions to / dpfvp/3, 

Jo 

where v is the cosmic-ray speed at momentum p. This gives 



n cr (r) 



6Pcr(r) ( m p c 2 



m n c 2 V 1 GeV 



B 1 a P — 2 3 - a } 



(88) 



where m p is the proton mass and B(z, w) = t z 1 {1 — t) w 1 is the beta function. Because 

Jo 

2 < a p < 3, the cosmic-ray pressure and energy density are finite even though the power-law 
scaling in equation (87) is taken to apply for p G (0, 00). Once n cr (r) is determined, the 



-25 - 



observed flux of gamma rays is obtained using equations (19) through (21) of Pfrommer & 
Ensslin (2004) and the model solutions for n e (r). The flux at Earth of photons with energies 
exceeding 100 MeV emitted by the central 300 kpc of each cluster is listed in table 2 for four 
values of the parameter a p , along with EGRET upper limits from Reimer et al (2003). The 
predicted fluxes are all below the EGRET upper limits. However, the fluxes equal or exceed 
the expected sensitivity to photons above 100 MeV of a future GLAST two-year all-sky 
survey [2 x 10~ 9 cm~ 2 s _1 ; Gehrels & Michelson (1999)]. Future GLAST observations will 
thus provide an important test of the AGN-driven convection model. 

Table 2: Predicted gamma-ray photon fluxes and EGRET upper limits 



Cluster 


Predicted flux above 100 MeV (10" 
a p = 2.1 a p = 2.3 a p = 2.5 


8 cm 2 s 1 ) 
a p = 2.7 


EGRET flux above 100 MeV 
(lO^cm^s- 1 ) 


Virgo 


0.36 0.55 0.45 


0.27 


< 2.18 


Perseus 


1.8 2.7 2.2 


1.3 


< 3.72 


A 478 


0.27 0.41 0.33 


0.20 


< 5.14 



Two sources of error in the predicted gamma-ray flux should be noted. As discussed 
previously, the model underestimates the plasma density within the central ~ 30 kpc, which 
acts to decrease the predicted gamma-ray flux. On the other hand, the inaccurate assumption 
that the cosmic rays are well mixed into the thermal plasma within the central ~ 30 kpc 
acts to increase the predicted gamma-ray flux. 

Hadronic interactions between cosmic-ray protons and thermal plasma also generate 
charged pions, which decay into secondary electrons and neutrinos. Pfrommer & Ensslin 
(2004) have calculated upper limits for the energy density of cosmic-ray protons in the 
Perseus cluster based on the radio emission from secondary electrons. Their upper limits are 
well below the peak values presented in figure 4, but depend upon the assumed magnetic field 
strength, the power-law index a p , and the spatial variation of p CT . This may be a problem 
for the model of this paper, and warrants further investigation. 
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4.6. Convection and the intracluster dynamo 

The origin of intracluster magnetic fields remains a mystery. If they originate from 
dynamo action in the intracluster medium, AGN-driven convection provides a means for 
generating the needed turbulence in cluster cores. 

I thank Eliot Quataert, Ian Parrish, and Jim Stone for discussions of the influence of 
parallel heat transport on convection, Steve Allen for providing the data for A478, Eugene 
Churazov for providing the data for the Perseus cluster, and Steve Cowley, Eric Blackman, 
Aristotle Socrates, and Olaf Reimer for other interesting discussions. This work was sup- 
ported by the National Aeronautics and Space Administration under Astrophysical-Theory- 
Program grant NNG05GH39G issued by the Office of Space Science. This work was also 
supported by DOE grant DE-FG02-01ER54658 at the University of Iowa. 
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